Efficient spin transitions in inelastic electron tunneling spectroscopy 
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The excitation of the spin degrees of freedom of an adsorbed atom by tunneling electrons is 
computed using a strong coupling theory. The excitation process is shown to be a sudden switch 
between the initial state determined by the environmental anisotropy to an intermediate state given 
by the coupling to the tunnelling electron. This explains the observed large inelastic currents. 
Application is presented for Fe and Mn adsorbates on CuN monolayers on Cu(lOO). First-principles 
calculations show the dominance of one coUisional channel, leading to a quantitative agreement with 
the experiment. 

PACS numbers: 68.37.Ef, 72.10.-d, 73.23.-b, 72.25.-b 



The way electrons flow through atomic contacts 
has important fundamental and technological implica- 
tions Electronic transport is a quantal process in 
which charge, spin and vibrational degrees of freedom 
are entangled leading to problems of intrinsic fundamen- 
tal interest. Technologically, the quest for minutariza- 
tion is pushing the limits of devices to the atomic scale, 
where the above transport properties will determine the 
actual device functionalities. An important issue is the 
appearance of inelastic effects where energy is taken from 
the electron flow into the different degrees of freedom of 
the system. Inelasticities lead to new regimes of trans- 
port that contain relevant information on the atomic con- 
tact and have been thus used to develop single atom and 
molecule spectroscopies 0, 0, Q • 

Inelastic electron tunneling spectroscopy (lETS) where 
electrons excite vibrations leading to conductance steps 
at certain voltage thresholds ^3 has been extensively 
studied in the last years d, [170, B Hi- The inelas- 
tic change in conductance is within a few percent of 
the elastic conductance, mainly due the smallness of the 
electron- vibration coupling 13, Recently, Heinrich 
and co-workers have been able to develop a spin-resolved 
spectroscopy using an STM 0, [l^, [l^, [3l • In magnetic 
lETS [1] , the tunneling electron yields energy to the spin 
of an adsorbed magnetic atom and in this way changes its 
orientation by overcoming the magnetic anisotropy bar- 
rier of the atom on the surface. Magnetic transitions in 
the meV range could be observed in adsorbates partly 
decoupled from a metal substrate p, l^, 13, UAUM- As 
in vibrational lETS, the conductance presents a step at 
the energy threshold however the changes in conductance 
at inelastic threshold can reach several hundreds percent. 
This is at odds with previous treaments [l3, H, 17| where 
first-order perturbation theory is used. 

In this letter, we present an all-order theory of the 
spin transitions lETS and apply it to the cases of Fe and 
Mn adsorbates on a CuN monolayer on Cu, experimen- 



tally studied in Refs. IJ, We compute the relative 
weights of both elastic and inelastic channels, leading to 
a quantitative account of the inelastic currents in the ex- 
perimental observations. The theory reveals the nature 
of the inelastic transitions and explains the extremely 
large inelastic currents in these magnetic systems. 

The general idea of our approach is the following. The 
spin of the adsorbate is in an initial state given by the 
anisotropy imposed by its environment and by an exter- 
nal magnetic field, B. During the very short coUisional 
time between the adsorbate and the tunneling electron, 
the electron spin couples with the adsorbate spin, forming 
a transient coUisional intermediate, whereas the interac- 
tion with the adsorbate environment can be neglected. 
This sudden switch between different coupling schemes 
of the adsorbate induces efficient transitions among mag- 
netic states. This excitation mechanism is not only found 
in STM-induced spin flip. Similar excitation processes 
have been shown to be very efficient for spin-forbidden 
electronic excitations in electron-molecule collisions 
or in surface processes [l^, as well as for rotational 
lETS 

The energy losses associated with the magnetic 
anisotropy in the presence of a magnetic field, B, have 
been modelled very efficiently in these systems 13|, |l4| 
using the following Hamiltonian : 



H = gfieB ■ S + DS'l + E{Sl - S^) 



(1) 



Where E and D are two constants describing the effect 
of the environment on the spin direction, g is the gyro- 
magnetic factor and /i b the Bohr magneton 1^ 14 1 . S is 
the spin operator of the adsorbate and Sx.y,z its projec- 
tions on the Cartesian axes. Diagonalisation of Hamil- 
tonian ^ yields the various possible states of the 
adsorbate spin in the system: 



\ct^n)=Y.'^r^,M\S,M) 
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FIG. 1: Electron transmission as a function of electron energy, 
tj, in a tunneling junction described by an atomic apex on a 
semi-infinite Cu(lOO) surface for the STM tip and a Fe (lower 
pannel) a Mn (upper pannel) atom on a CuN monolayer on 
a semi-infinite Cu(lOO) electrode. The magnetic atom-apex 
distance is 5.2 A in the present calculation. The full line is 
the transmission for the majority spin and the dashed line for 
the minority one. 



where \ S, M) are eigenvectors of the S*^, Sz operators. An 
electron injected from the STM tip collides with the ad- 
sorbate on the surface and can cause inelastic transitions 
between the (j)n states, which are recorded in an lETS 
experiment. In the present work, we use the modelling 
performed in Ref. [ij, (D and E parameters, g and 
spin of the adsorbate), which very precisely reproduces 
the energy positions of the inelastic thresholds. 

In a first step, we compute the electron transmission 
through Fe and Mn adatoms on a CuN monolayer on 
Cu(lOO) by DensityFunctional Theory (DFT) with the 
Transiesta code [2l|. We used the generalized gradi- 
ent approximation (PBE) and a double-zeta local basis 
set where the contact region is modeled by a 7-atom 
slab, a CuN layer and a Fe (Mn) atom, a vacuum gap 
of 5.2 A and a 5-atom slab with an extra atom for the 
tip region of the contact. The contact is relaxed using 
the Siesta method 
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^ ^ Atomic forces are relaxed be- 
low 0.04 eV/A. The transmission is then computed for 
zero bias voltage, using the bulk Cu unit cell along the 
[100] direction as the primary unit of the two semi-infinite 
electrodes coupled to the contact region. On the energy 
scale relevant for the present magnetic lETS context, the 
majority spin transmission is 20 times the minority spin 
one for the Fe junction, see Fig. [1] Despite the fact that 
the actual spin state of the atom cannot be taken into 
account by DFT, these simulations can yield quantita- 
tive data in spin transport [l^. Indeed, this difference 
implies that the transmission of minority-spin channels 
can be considered as suppressed. In addition, the trans- 



mission is seen to be almost flat as a function of electron 
energy, so that simple branching ratios can be used to 
obtain the relative value of the elastic and inelastic con- 
ductance of the system as a function of the STM bias. For 
the Mn adsorbate, we find a similar result, though with 
a weaker dominance (factor 5) of the majority channel 
(Fig.©. 

The branching ratios between elastic and inelastic con- 
ductance are determined making use of the following 
facts: 

i) the rotation of the adsorbate spin, S*, due to the 
magnetic anisotropy, Eq. ([T]), is slow compared to the 
electron-atom collision time so that we can use a sudden 
approximation, neglecting the effect of Hamiltonian ([1]) 
during the collision. 

ii) the spin of the tunnelling electron couples to the 
spin of the atom to define collision channels of total spin 
5t = S" + 1/2 and S* — 1/2 that are linked to the asymp- 
totic channels of the collision via: 

\St, Mt) = CGs^,MT,m \S, M = Mt - m) |l/2, m) 

m 

(3) 

where the kcts on the rhs correspond to the decoupled 
spins of the atom and of the timnclling electron, m is the 
projection of the electron spin on the z-axis. The CG 
are Clebsch-Gordan coefficients that give the weight of 
the various \S,M) states in the collision channels. From 
Eqs. ([2]) and ([3]), we can express the collision channel 
states as functions of the initial and final states of the 
collision: 

|j) = |5t, Mt) = J2 I0n> |l/2, m) (4) 



It yields the weight of the various anisotropy states in 
the collision channels associated to the total spin, St- 

iii) From the DFT result, we only consider the maxi- 
mum spin intermediate state {St ~ 5/2) for the Fe ad- 
sorbate and {St = 3) for Mn adsorbate. 

iv) From Eqs. ([4]), we can derive the amplitude for 
transitions from |0„)|l/2,m) to |0„') |l/2, m') through 
the intermediate j as proportional to the product 
Aj^n,m^j,n',m' ■ The Contributions from the different in- 
termediate states are then added coherently for the in- 
distinguishable channels (same final (n',m') state for a 
given (n, m) initial state) and incoherently for the distin- 
guishable channels leading to the relative excitation prob- 
ability (branching ratio) of the different excited states: 
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FIG. 2: Computed conductance for a Fe atom on a CuN 
monolayer on Cu(lOO) in atomic units. The conductances for 
increasing magnetic field B = are vertically displaced for 
representation purposes. The B field is oriented along the N 
axis in part (a) and along the hollow axis on the surface in 
part (b). 



Note that because of the dominance of one St interme- 
diate state in the conductance, the sum over j only con- 
cerns the M sublevels, i.e. the orientation of the spin of 
the intermediate state; the corresponding contributions 
only differ by spin coupling coefficients and add coher- 
ently. The sum over j runs over the St = 5/2 (resp. 
St ~ 3) intermediates for the Fe (respectively Mn) ad- 
sorbates, and the sum over m and m' concerns the spin 
up and down of the coUisional electron. Equation ([5]) 
above has been derived for an unpolarised incident elec- 
tron; it can be easily generalized to yield spin-resolved 
transitions. 

Equation ([5]) is the basis of the present work. It yields 
the relative weight of the clastic and inelastic channels in 
the conductance. This expression is a direct consequence 
of spin coupling and magnetic anisotropy, associated to 
the dominance of the majority spin conductance. 

Fig. [2] presents the conductance as a function of the 
STM bias obtained as the product of the computed global 
conductance (Fig. [T]) by the elastic/inelastic branching 
ratio from expression Results arc shown for Fe ad- 
sorbates at five values of the B field (B along the N axis 
in part (a) and along the hollow axis in part (b)). A 
gaussian broadening of 0.26 meV corresponding to a tem- 
perature of 0.5 K 2J| has been added. In this system, 
the Fe spin is equal to 2 [l^l so that the conductance 
can present 4 steps associated to the inelastic thresholds 
(labelled 1-4 on the figure). As a first remark, the con- 
tribution of inelastic channels is very large; for B = 
and for an infinite resolution in this system, the inelastic 
channels at large bias amount to around 67% of the elas- 
tic channel. At finite resolution, for B=0, the increase of 
the conductance between and 8 mV is smaller due to 
the small energy difference between (j)Q and 0i. Second, 
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FIG. 3: (Color online) Relative inelastic step heights in the 
conductance for Mn and Fe adsorbates as a function of the 
applied magnetic field, B, (along the N axis): calculations 
(full lines) and experiment (symbols). 



Fig.[2]shows an important change in the inelasticity spec- 
trum with B. The 0-1 excitation is dominating at low B 
and disappears when B increases, whereas 0-2 dominates 
at large B. The 0-4 excitation is always weak. This be- 
haviour is exactly the one observed experimentally [ist . 
For a quantitative comparison. Fig. [3] presents the rela- 
tive step heights (ratio of the height of a given inelastic 
step to the sum of the inelastic steps 1-3) as a function of 
B, compared with the experimental values. The 0-4 exci- 
tation is predicted to be very small and it is not observed 
experimentally for this geometry; we have not included it 
on the figure. Results obtained for the other orientations 
of the B-field also reproduce the importance of inelastic 
channels. 

Experiments on Mn adsorbates on CuN monolayers on 
Cu |13| showed a very small magnetic anisotropy asso- 
ciated with a spin 5/2 and at finite resolution, the con- 
ductance is basically exhibiting a single inelastic step for 
all B values. Fig. [3] presents a comparison of our pre- 
diction for the relative inelastic step height (ratio of the 
inelastic step to the conductance at bias) as a function 
of B, it is seen to be in quantitative agreement with the 
experimental data. 

In the present approach, the excitation process is seen 
as a decoupling/recoupling process induced by the col- 
lision with the tunnelling electron. As our theory pre- 
dicts, spin excitation can take place without spin-flip of 
the electron flux. We can distinguish two regimes as the 
magnetic field, B, increases. These are associated with a 
change in the magnetic structure of the system: evolution 
from a magnetic anisotropy induced by the lattice at low 
B, towards the Zceman effect at higher B. The main axis 
of the lattice- induced anisotropy is the N axis (z-axis) and 
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thus, the efficiency of the B field in generating a Zeeman 
structure is weaker along the hollow axis, as compared to 
the N axis (see Fig. [5]). This change of the energy land- 
scape has bearings on the actual spin composition of the 
atom states. As an example, the 0-1 excitation of Fe, see 
Fig. [21 comes from the decoupling with the environment 
at low B, and at large B, the excitation is mainly due 
to coupling to the electron spin. Furthermore, in Fig. [3l 
for incident electrons polarized along the N axis, the 0-1 
and 0-4 transitions are not associated to a change of the 
coUisional electron spin direction, whereas the 0-2 and 
0-3 transitions are entirely spin-flip transitions. 

In the case of Mn, the environment-induced anisotropy 
is very weak and for finite B, the Mn spin structure is 
a simple Zeeman splitting; in this case, the transitions 
are only spin- flip with a AM = ±1 selection rule and 
the fraction of inelastic tunnelling is basically given by a 
ratio of squared Clebsch-Gordan coefficients leading to a 
nice agreement with experiment, see Fig.[3l 

The importance of the inelastic conductance in the to- 
tal conductance is a direct consequence of the nature of 
the excitation process, analysed above. The initial mag- 
netic state of the adsorbate couples with the spin of the 
collisional electron to form a collisional intermediate with 
a given total spin. At the end of the collision, the colli- 
sional intermediate populates all the possible asymptotic 
channels according to their weight in the intermediate, 
Eq. ([5]). The importance of a particular inelastic channel 
is then given by its weight in the collision intermediate 
(environment-induced anisotropy or spin coupling coeffi- 
cient) and is not proportional to the modulus square of 
a matrix clement between initial and final states. This 
explains the large difference of inelasticity observed in 
magnetic lETS compared to vibrational lETS. 

However, inelastic rotational excitation shares many 
features with magnetic lETS. The resonant rotational 
excitation of an adsorbed molecule induced by tun- 
nelling electrons also involves a transient angular mo- 
mentum coupling between the collisional electron and the 
molecule. The corresponding process can be formulated 
in a very similar way to the present work and indeed leads 
to the observed strong rotational excitation 2^, 2^. 

We believe that the present mechanism, which explains 
the strength of magnetic lETS and accounts for the ob- 
servations in the case of Fc and Mn adsorbates, is of 
general occurrence. In addition, the present formalism 
yields a very easy way of accuratly predicting the im- 
portance of spin transitions. One can stress that, in the 
present approach, once the spin state of the collisional in- 
termediate is fixed, all the spin changing transitions are 
fully determined. Indeed, the above formalism also leads 
to the quantitative account [i^l of the magnetic lETS of 
Mn2, Mus and Co-phthalocyanine layers In the 
latter case, we find an extremely strong inelastic frac- 
tion in tunnelling, up to 300% for 3 molecular layers, in 



excellent agreement with the experimental findings. 
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